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1. Introduction 



In the Lattice2007 symposium we showed [|l|] that the infrared QCD running coupling of 
ghost-gluon vertices of lattice Coulomb gauge fits better than those of lattice Landau gauge with 
the experimental running coupling data extracted from the quark-gluon vertices The gauge con- 
figurations of Kogut-Susskind (KS) fermions of the MILC collaboration [|J] and those of Domain 
Wall Fermion (DWF) of RBC/UKQCD collaboration [Q] were used in these simulations. 

Recently I found that the quark propagator of DWF can be calculated by an extension of the 
conjugate gradient method that we adopted in the calculation of the quark propagator of KS fermion 
or Wilson fermion [§]. I calculate the Coulomb gauge quark propagator near the cylinder cut and 
small momentum transfer using the DWF gauge configurations and extract the quark-gluon vertex, 
which is written as 



r n(p,l) = S 1 (p)Gff(p,q)S \p,q) 



S ab [gi(p,q)Yn + ig2(p,q)qfi +g3{p,q)p^4 



(1.1) 



as scalar, vector and tensor form factor, respectively. A general method for non-perturbative renor- 
malization of lattice operators was presented by Martinelli et al.0 and I adopt their method to 
extract the Coulomb gauge form factors. The lattice simulation of DWF is given in [JsJ]. In the 
mass-less limit, the wave function possesses exact chiral symmetry and its non-perturbative renor- 
malization is investigated [|, 1C]. 



2. The lattice Coulomb gauge 



The minimizing function of the Coulomb gauge(d,A,- = 0) in the Log-U version [ 13] is Fy [g] = 
I \A S 1 1 2 = L A y tt {^x/Ax i) ■ R emnant gauge fixing is not done although the gauge field Aq(x) can 
be further fixed by the following minimizing function via the gauge transformation g(xo). 

The gauge configurations that I adopted are summarized in Table 1 . In the process of conjugate 



N/ 



m 



l/a(GeV) L, L, aL s (fm) 



DWFoi 2.1308/) 2+1 0.01/0.04 1.743(20) 16 32 1.81 
DWF 02 2.1303/) 2+1 0.02/0.04 1.703(16) 16 32 1.85 
DWF 03 2.13(0/) 2+1 0.03/0.04 1.662(20) 16 32 1.90 



Table 1: The parameters of the lattice configurations. 



gradient iteration, I search the shift parameter for (j>i and for 0/? as follows. In the first 50 
steps I choose au = Min(afc, af) and shift = (p^ — ttkQfc and <pj* +l = tyf — d^tyf and in the last 
25 steps I choose at = Max(aj:, a*), so that the stable solution is selected for both (j) L and (f> R . 

The convergence condition attained in this method is about 0.5 x 10~ 4 . One can improve the 
condition by increasing the number of iteration, but the overlap of the solution and the plane wave 
do not change significantly. 

In our Lagrangian there is a freedom of choosing global chiral angle in the 5th direction, 

¥ _» e m Wi y e ~ ims y. (2. 1) 
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m u d/a m s /a c A(GeV) a 



DWFoi 


O.Ol 


0.04 


0.49 


0.76(2) 


1.25 


DWF 02 


0.02 


0.04 


0.48 


0.80(3) 


1.25 


DWF 03 


0.03 


0.04 


0.61 


0.66(2) 


1.25 


MILC/i 


0.006 


0.031 


0.45 


0.82(2) 


1.00 


MILC /2 


0.012 


0.031 


0.43 


0.89(2) 


1.00 



Table 2: The fitted parameters of mass function of DWF(RBC/UKQCD) and KS fermion (MILC). 

On the lattice, the expectation value of the quark propagator S(p) consists of spin dependent 
srf p part and spin independent 88 part. I specify the propagator of the left-handed quark by the 
suffix L and the right-handed quark by the suffix R. With use of the renormalization factor Z, they 
are defined as 



and 



Tr(x(p,s)P L/R ^(p,s))=Z(p)(2N c )m L/R (p,s 



Tr(x(p,s)ipP L/R *¥{p,s)) = Z(p)/(2N c )ip*f L/R (p,s), 



where p\ = - sin — — (p; = 0, 1,2, ■ • • ,«,/4). The fifth coordinate of the DWF is specified by 
a rii 

5 = 0,1, •••,15. 

We parametrize the mass function j$ = 88 '/ ' stf as 

I tried a = 1, 1.25 and 1.5, compared % 2 and found that a = 1.25 gives the best global fit. Since 
the pole mass Q W is not included in these plots, m / is set to be here. 

In the case of KS fermion of m f = 0.0136GeV and 0.027GeV data [@ ||], we fixed a = 2 



and obtained A = 0.82GeV and 0.89GeV, respectively. In general A becomes larger for larger a, 
but A of DWF seems larger than that of KS fermion. In the case of KS fermion, A becomes smaller 
for smaller mass mf, but in the case of DWF, it is opposite. 



3. The vector current quark-gluon vertex 

Near the cylinder cut and when the momentum transfer q is small, the quark gluon vertex in 
momentum space is calculated from 




where 5; is a DWF propagator of the i'th sample among altogether Af samples. 

I specify the spin part of a/3 as H — which corresponds to spin up for a and spin down for 
j3. In the case of T = q\Jiqu I choose 74/i + ~ or 7472 H as an example and £/j(p\0)~ + as the 
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p(GeV) 

Figure 1: The mass function of the domain wall 
fermion as a function of the modulus of Euclidean 
four momentum p. DWFoi. (149 samples). Blue 
disks are ml (left handed quark) and red boxes are 
niR (right handed quark). 



53 



q[GeV] 

Figure 2: The running coupling 0£ ligl (q) of 
MILC/i (blue disks) and DWFoi (green points). 
The dash-dotted line is the pQCD result and the 
dashed line is the pQCD with the (A 2 ) condensate 
contribution. The red points are the experimental 
data of the JLab group. 



incoming wave £/i(p\0) ++ " f contributes as the outgoing wave. When jz^(/?|0) is the incoming 
wave £^(p|0)~ + ^ contributes as the outgoing wave.A@ 
The quark propagator is given as 



S(P) 



-istf p + , 



(3.2) 



I evaluate vector current matrix elements by diagonalizing the matrices in the eq.(3.1) and getting 



the j2 of the trace. 



3.1 The scalar form factor and the running coupling 



The vector current Ward identity allows us to extract the running coupling a Sigl (q) from the 
difference of (p + f ) and S~ l (p - § ) [|] 



-i[S-\{p + \)^)-S-\(p- q -) j \0)]=Z v K (p) qj /AK. 



(3.3) 



When a crossing is performed and a quark is transformed to anti-quark, p can be treated as the 
momentum transfer q, to be compared with the experiment. The green points above q = 1 .4 GeV 
in Fig.^| are the result of DWF. The error bars are taken from the Bootstrap method after 5000 re- 
samplings [0, [15] ], which are smaller by about a factor of 10 as compared to the standard deviation 
of the bare samples. The data are comparable with the pQCD results with a phenomenological 
(A 2 ) condensates effect observed in the ghost-gluon coupling of the MILC/i Landau gauge fixed 
configuration [13]. Although statistics is not large, I observe that the running coupling a sgl (q) 



of Landau gauge fixed configuration of MILC/i is not infrared suppressed in contrast to a s (q) 
measured from the ghost-gluon vertex [Oj. 

I compare the lattice results with experimental results extracted by the JLab group [||, [||. The 
JLab group analyzed the difference of the spin-dependent proton structure function and the spin- 

Q 2 



dependent neutron structure function as a function of x 



2Mv 



where Q is the four-momentum 



squared, V is the energy transfer and M is the nucleon mass. The difference reduces the contri- 



bution of A(1232) and at infinite momentum transfer squared the Bjorken sum rule [16], and at 
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Figure 3: The tensor form factor h.2{p)p4q/ pn\q\ 
at \q\ = 2.6GeV. (DWFqi, 52 samples) 
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Figure 4: The vector form factor A\(p)q/\q\ at 
\q\ = 2.6GeV. (DWFqi, 52 samples) 



momentum transfer the Gerasimov-Drell-Hearn sum rule [17, 18] allows to fix the renormalization 



[19]. They extracted the running coupling in the infrared from 

a, 



r? 



f[l-^ + (a>)]. 

6 K 



(3.4) 



In a confining theory where gluon have an effective mass, all vacuum polarization corrections to 
the gluon self-energy decouple at long wave length and one expects an infrared fixed point. A 
closeness of F? and F" implies a presence of infrared fixed point a s gl (0) = K. When the coupling 
is constant in the infrared and the quark masses can be ignored, one expects conformal symmetry 
or ADS/CFT correspondence to be applied in non-perturbative region j pfj| , [2l| ]. Analytical Dyson- 
Schwinger approach also suggests the presence of infrared fixed point [g^]. 

3.2 The tensor form factor 

I specify the spin of quarks that sandwich 7471 , 7472 and 7473 and pick up the quark propagator 
Si(p\0) whose spin of the final state matches with that of 747. I assume propagators are color 
diagonal. The vertices 7471 and 7472 are spin off-diagonal and the vertex 7473 is spin diagonal. 
There is a relative phase factor i in 7472 as compared to the other 747. 

The tensor term is evaluated from the difference of 



«(P ~ |)74 ~))P4(P+~)j ~ (■^(p)7A J £7j^/(p))P4Pj 

I sample-wise diagonalize 

F L / R = [{^ L l\(p\0)p4(Y 5 ^ i (pW^Y 5 )ofp l )o l 

+ (^ L/R 4(P\0)P4(75^ L/R i(p\0) a ^r5)0^p2)0 2 

H^ LIR MQ)(j^ LIR i{pW^Y5)°?p^] 



(3.5) 



(3.6) 



for each a,j8 that specify the color and the spin, and take the sum of the real part of the positive 
eigenvalues. The corresponding momentum shifted matrix elements are 



pL/R 

A A 



4 ,R {P-\\^)P^ 



t'\p + \W^K)°f{P + \)^ 



+K i/ ^P-||0)^4(75^b + ||0)^ t 75)a 2 a/3 b + |) 2 )a2 
+K L/R b-||0)^(7 5 ^ L/R b + ||0) a/3t 7 5 )a 3 a/3 b + |)3)a3]. 
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I approximate the denominator of the propagator for Fa and by those of the quark propagator in 
the cylinder cut and evaluate the tensor term as 



i (fj+f*)-(r/+iV^ 

12 {sip>- + Jl3g){si*iP- + JtL 



l K±AZ±A> \LA Z±A I = z v K 2 (p)p, qi /AK (3.7) 



A simulation result of DWFoi is shown in Fig.[|. 

3.3 The vector form factor 

I define the matrix elements between spin ++ as Yc and between spin as T^. 

r c ,j = Re[(^ L++ {p-^\0)y A {y 5 ^ ++ { P + q -j-WY 5 ) 
_^++ (p _||0) 74 ( 75 ^++(p + ||0)t 75 ))] 

and 

T DJ = Re[(^--{p-^\0)y A y 5 ^--{ P+ q -l\0^y 5 ) 

_^-( p _||0) 74 ( 75 ^-(^ + ^|0) t 7 5)}] (3-8) 

- E \ T D,j\ + \Tc,j\ = z y {) /AK 

A simulation result of the vector form factor of DWFoi is shown in Fig.|j. 
4. Discussion 

I showed that nonperturbative renormalization of lattice Coulomb gauge quark-gluon vertices 
is feasible. In addition to the running coupling, I measured the tensor term g3 (/>,<?) 74/?4^ and the 
vector term g2(p,q)q of the Coulomb gauge quark-gluon vertex. 

The running coupling (Xi(q) and a s , gl (q) of DWF are consistent with the JLab extraction. The 
running coupling of Landau gauge a^, (q) is not infrared suppressed in contrast to a s {q) of triple 



gluon vertex [ ]23| ] and ghost-glion vertex [JTJ, |TJ]. The infrared suppression of tt s {q) in Landau gauge 
was attributed to an instanton effect [p3|]. However, as shown by 't Hooft [p4]], in the instanton 
calculus extended to the infrared region, there is a divergence from zero modes. In supersymmetric 



case, however, the divergence from fermionic zero mode and bosonic zero mode cancel [ 25 ] and 
a finite infrared fixed point appears. It is plausible that the search of a stable solution in Coulomb 
gauge with use of the conjugate gradient method leads to a solution in which zero mode divergences 
of gluons and those of quarks cancel. 

When DWF configurations of larger lattices are available, the a s , gl (q) data can be extended to 
lower energies. I expect that CC Sigl (q) of DWF in Coulomb gauge show the same behavior as CCj(q). 
The difference of the infrared features of Landau gauge and Coulomb gauge running coupling casts 



questions on the Kugo-Ojima confinement criterion [126, 27] which is derived in Landau gauge but 
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effects of instantons are not considered. Simulations with momenta far from the cylinder cut and a 
calculation of the form factor of three quark systems are also left in the future. 
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